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Abstract 



We analyze the long time behavior of transport equations for a class of dissi- 
pative quantum systems with Fokker-planck type scattering operator, subject to 
confining potentials of harmonic oscillator type. We establish the conditions 
under which there exists a thermal equilibrium state and prove exponential de- 
cay towards it, using (classical) entropy-methods. Additionally, we give precise 
dispersion estimates in the cases were no equilibrium state exists. 
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1 Introduction 

In this paper we analyze a class of dissipative quantum systems, modeling the 
motion of a particle ensemble, say electrons, interacting with a heat bath of 
oscillators. The resulting irreversible dynamics for the electrons, is a typical 



example of a so-called open quantum system (Da|, i.e. a system in which the 
interaction with the environment is taken into account. The evolution equation, 
sometimes called Master equation, for the density (matrix) operator R(t) of the 
particles reads 

±R = -~[H,R]+A(R), (1.1) 
R(t = Q) = R , (1.2) 

where H is the self-adjoint Hamiltonian operator of the free system and A 



models the effects introduced by the heat bath. The right hand side of fll.l| ) 
constitutes the formal generator L of the quantum dynamical semigroup acting 
on R. 
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Assuming that the time evolution of the system is Markovian, G. Lindblad B 
gave the most general form of a bounded operator L, such that the semigroup 
preserves the positivity, hermiticity and the normalization (unit trace) of the 
density operator R. However, for unbounded operators L, which is the case 
in our work, the so-called Lindblad condition is necessa ry b ut not sufficient to 
guarantee the conservation of these properties (see, e.g., [CF] and the references 
therein). 



Using the Wigner transform [Wi], dissipative quantum models can be equiva- 
lently represented in phase space, resulting in a kinetic transport equation with 
interaction terms for the quasiprobability distribution of the particles. In this 
paper we assume that the mechanism coupling particles and environment can 
be described by a linear scattering operator C q of Fokker-Planck type. The 
Lindblad condition therefore reduces to the assumption of a positive definite 
diffusion matrix D. 

Up to now, a mathematically rigorous derivation of such a Quantum Fokker- 
Planck equation (QFP) from many-body quantum mechanics is still missing. To 
the authors knowledge, the only result in this direction is given in [CEFM], which 
however justifies only a particular case of the class of models considered in this 
work. Nevertheless there exists a huge amount of a somewhat phenomenological 
physical literature on this type of equations, which play a relevant role within 
the areas of quantum optics (laser physics) ]Da|], [ Dc 
[Ht|, quantum brownian motion ]CaLe |, [ pi[ | , |HuMa 



microelectronics 
and the description of 



decoherence and diffusion of quantum states [AnHa|, |DGHP|. 

Rigorous well-posedness and existence of local in time solutions of the frictionless 
QFP equation, with self-consistent Coulomb interaction, have been studied in 



a precedent paper of one of the authors [ALMS|. The present work investigates 
the long time behavior of the linear QFP equation in the presence of friction and 
an exterior time-independent potential V. The phase space formalism provided 
by the Wigner transform proves to be particularly helpful for this task, since it 
allows the use of certain entropy techniques established for cla ss ical diss ipative 
equations (for an overview on these techniques, see, e.g., | MaVi ] , AMTU |). The 
word "entropy" is used here in a mathematically technical sense and can be seen 
as a generalization of the classical entropy concept of L. Boltzmann. It should 
not be confused with the physical correct von Neumann entropy of quantum 
states. 

By comparison with the classical Fokker-Planck equation (FP), we expect the 
solution of the QFP equation to approach a thermal equilibrium state in the 
long time limit, provided the friction term appearing in C q is positive and the 
exterior potential V is confining, i.e. V(x) — > oo as |a;| — > oo, fast enough. 
In our work we shall specify the potential V to be harmonic. This particular 
choice allows explicit calculations and is maybe the most fundamental one, from 
a physical point of view [ CEFM ] . Using the entropy approach, we will prove the 
convergence of the solution towards the steady state with a precise exponential 
rate, under the assumption that the initial data wq has bounded "entropy," 
relatively to the equilibrium state. 

This paper is organized as follows. In section 2 we set up the model and collect 
some preliminaries. In section 3 we specify the potential V to be of harmonic 
oscillator type and explicitly calculate the corresponding equilibrium state. Ex- 
ponential decay towards it will be proved in section 4, where we also give precise 
dispersion estimates in the unconfined cases. 
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2 The model: preliminaries 

We consider a linear dissipative equation modeling the motion of particles, say 
electrons, under the influence of an electric scalar potential V and a thermal 
bath of harmonic oscillators in thermal equilibrium. In the following we denote 
by pit,-) £ L 2 (M. 2d ) the N-particle density matrix function of the electrons, 
which is the kernel of the self-adjoint trace-class density (matrix) operator R(t) : 
L 2 (R d ) -> L 2 (R d ), i.e. 



(R(t)f)(x) := / p(x,y,t)f(y)dy. (2.1) 
The evolution equation of the electrons is given as a PDE for the density matrix 



d t p 



-{H x -H y )p~ 1 {x-y) ■ (V a 



v) 



D qq \V x 



V y \ 2 -^\ X -y\ 2 + 2 -^L 



subject to the initial condition 

p(t = 0,x,y) = po(x,y), x,yeR d . 
Here H x (rcsp. H y ) denotes the electron Hamiltonian 



H ,■ 



— z 



V{x), 



{x-y)-(V x + V v )\p, 
(2.2) 

(2.3) 
(2.4) 



acting on the x (resp. y) variable. The constant m stands for the mass of the 
individual particles. Equation (2.2) is a generalization of the Caldeira-Leggett 
master equation for medium temperatures, see f5f |, [Dc]. 

On a kinetic level this model reads (QFP equation) 



d t w + £ ■ W x w + e[V]w = C q 



w(t = 0,x,£) = w (x,£), 
where the scattering operator C q is defined by 

£ w ~ ^-Atw + 27dW£u>) + D qq A x w + 2^ div x (V £ w). 
m z m 



(2.5) 
(2.6) 



(2.7) 



Here w(t,-) £ L 2 (M. d x M$) is the Wigner transform Wi] of the corresponding 



density matrix pit, ■) S L 2 (R d x 
1 



w(x,£,t) 



{2n) d h 



i> ( x + - — y, x — - — y, t ) e l ^' v dy. 
2m 2m 



(2.8) 



In the literature w(t, •) is often referred to as a quasiprobability distribution 
because it generally assumes negative values too [ |pq| , [Hu|. 
The (non-local) pseudo-differential operator <d[V] is defined by 



Q[V]f(x,0 := 



V[x+—y 

2m 



V x 



y 



2m 

f(x,C) e^-^d^dy 



(2.9) 
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In order to be consistent with the usual density matrix formulation of open 
quantum systems |Da] in the class of Lindblad operators O], we assume for the 
diffusion constants D pp > and D qq ,D pq > 0. On the other hand, the fric- 
tion parameter 7 has to be nonnegative. Additionally we impose the following 
relation (for more details see [ALMS], |Q) 



DppDqq 



D 2 pq >— - and £ pp >0if 7 



0. 



(Al) 



Note that if 7 > 0, condition (Al), the so-called Lindblad condition, implies 
that (2.7) is a uniformly elliptic operator. 

Using the Wigner transform the charge and flux densities associated to the 
density matrix p can be defined (formally, since w $ L 1 in general) in the same 
way as in classical statistical mechanics. Namely they are given as moments of 
the Wigner transfom, (for details see, e.g., [GaMa], [LiPa]), 

n(x,t) := / w(x, £, t)d£, 

j{x,t) := / £w{x,£,t)d£. 



(2.10) 
(2.11) 



Although not obvious from the above definition, the necessary positivity of n is 
guaranteed for physical quasiprobabilities w, i.e. for w's, which are indeed the 
Wigner transformed kernels p of density operators R (see [Ar , [LiPa], [Ta] for 
a more complete account on this). 

With the above definitions we obtain, after formally integrating the QFP equa- 
tion (2.5) w.r.t. £, the associated "continuity equation''' 

d t n + div j = D qq A x n, (2.12) 

which obviously implies (again on a formal level) the conservation of mass, i.e. 



M(t) := 



w(x, £, t)dxd^ 



In view of this property, we assume for simplicity that 



M n = 



Wo(x,!;)dxd£ = 1. 



(2.13) 



(A2) 



This is of course not a restriction as long as the equation is linear. The continuity 
equation (2.12) suggests that in this dissipative model the usual definition of 
the flux density (2.11) needs to be replaced by 



J(x,t) : = j(x,t) - D qq V x n(x,t) 

= / {i~D qq V x ) w{x,S„t)dt 

JR d 



Thus, instead of ( 2.12 ), we obtain the following conservation law 
dtn + div J = 



(2.14) 
(2.15) 

(2.16) 



associated to the QFP equation. 
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Remark. In physical units, the friction and diffusion constants are usually 

XH 2 



given by (see, e.g., pD, Q) 



' 2m' pp ' 99 12m 2 fc B T' 



P9 ~~ 127rmfc B T' 



(2.17) 



Here A > is the coupling constant of the heat bath, fc B the Boltzmann constant, 
T the temperature of the bath and Q the cut-off frequency of the reservoir 
oscillators. In terms of these constants the Lindblad condition (Al) reads 

k B T 



n < 



(2.18) 



which implies the validity of our model at medium/high temperatures. 

In the classical limit h~ — > 0: D qq = D pq — and (at least formally) the pseudo- 
differential operator simplifies to 

1 



®[V]f 



-v x v ■ V 6 f. 



(2.19) 



Thus we recover the well known kinetic Fokker-Planck equation for the limiting 



classical phase space probability distribution w {t, •) G 
of positive bounded Borel measures) 

1 



i) (the cone 



d t w cl 



-V X V ■ W 5 w cl = 



'-A 6 w cl + 2 7 div<;(£w cZ ). (2-20) 



Also, the flux density J formally simplifies to the classical one, i.e. J — > j, as 
h — > 0. For a rigorous theory of such homogenization limits see, e.g., |LiPa], 
[GMMP] and for details on the extensively studied FP equation, see for exam- 
pie QRi|. 

Furthermore note that the Lindblad condition (Al) disqualifies the classical 
FP scattering operator [CEFM| as a relevant quantum mechanical model for 
the environment interaction. 

In this work we are concerned with the following solution concept for our IVP. 



Definition 2.1. A function w <E C(R^;L p (Rf x Rf)), with 1 < p < oo, is a 



mild solution of the IVP (fi.dj , \2.(\ ) with w a G i p ( 



1* x Mi), if and only if 



■](t,x,£) = 



R d xR d 



wo(xo,£o) G(t,x,£,x ,£o) dx d£a, 



(2.21) 



where the Green's function G satisfies equation (2.5) for all fixed (xoj£o) ^ ^ 2d , 
all t > 0, with an initial condition 



\imG(t,x,£,x ,£o) = 5(x- xo,€-€o). 



(2.22) 



that has to be understood in a weak sense. If additionally w G C 1 (R + ; C 2 (M^ x 
Rf) 



p< ?)), then w is called a classical solution. 



Remark. In classical physics the theory of kinetic equations focuses on L 1 - 
solutions. This implies mass conservation, since classical phase-space distribu- 
tions are pointwise positive functions. Having in mind the definition of the 
Wigner transform, the L 2 -norm is more convenient in our quantum mechanical 
context. 
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3 Harmonic oscillator potentials 

In the next two sections we choose a normalization such that K = m = 1, for 
simplicity. We moreover assume the confining potential to be of the following 
class 

uj 2 

V(x) = -^-\x\ 2 + ax + b, a,beR,ujQ>0. (A3) 

An easy calculation shows that, maybe after an appropriate shift in the in- 
variable, the pseudo-differential operator Q[V] is given by 

Q[V]w = -to 2 x ■ V ( w. (3.1) 

The QFP equation thus simplifies to 

d t w + £ • V x w - w%x ■ V s w = £ g w x,£ e M d ,t e M.+ (3.2) 
w(t = 0, x, £) = w (x, £), (3.3) 

which can be equivalently written in the more compact form 

d t w = div( XiS ) (DV( x tfW + P(x, g)w) , (3.4) 
w(0,x,£) = w (x,0, (3.5) 

where the diffusion matrix D and the vector-valued drift P are given by 

D qq i d D pq i d \ Df _ rt _ ( -e 



Here Id denotes the idendity matrix in R . Note that the Lindblad condition 

ite 

2d 



(Al) guarantees that the diffusion matrix is positive definite if 7 > and then 
( O) is parabolic in the phase-space coordinates (x, £) S 



3.1 Fundamental solution 



Now consider only the first order part of the operator (3.2), resp. (3.4). The 
associated characteristic ODE's are given by 

X = £, X(t = 0)=x , (3.7) 

e = -H 2 X + 2 7 £), £(i = 0)=$>. (3-8) 

This system defines the characteristic flow $4(2:0, £0) = [Xt(xo, £0), X t (xo, £0)] 
in phase space if x 1^. This flow can be explicitly calculated, depending on 
the size of the friction constant 7. 



Lemma 3.1. Consider the system (3.7), (3.t ) with 7 > 0, then the dissipative 
flow $ t : R 2d — > R 2d is given by: 



1- If < 7 < ujq, then, with to := \J uJ^ — j 2 , 

e -7* 

$t(zo,£o) = [x (ucos(ut)+jsm(ujt))+£ Q sm(ujt), 



£o(w cos(wi) — 7 sin(cjt)) — xqWq sin(wi)] . (3.9) 
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2. Ifj>ujQ, then, with u> := y^j 2 — lu 2 , 



\xo(u> cosh(wt) + 7sinh(W)) + £o sinh(wt), 
£o(u> cosh(wi) — 7 sinh(ijji)) — xqujq sinh(wi)] . 



(3.10) 



3. If 7 = loq then 

$t(a?o, &) = [( 7 t + l)zo + f£o, (1 - lt)to - l 2 tx ] . (3.11) 

Proof. The proof follows from straightforward calculations. □ 

Now, using lemma 2.1, we obtain an explicit representation of Green's func- 
tion G of the QFP equation with harmonic oscillator potential, depending on 
the size of 7. 



Proposition 3.1. Let 7 > and let conditions (Al), (Ah) hold, then the 
Green's function G associated to (S.i) is, for every fixed t > 0, a pointwise 
positive function, given by 



G(t, x, £, x ,to) = e 2d ~< l F{t, X_ t (x, - x 0) X_ t (x, - 



(3.12) 



with 



l ' ' U (2 7 r)' i (4A(<)^(t)-/i 2 (t)) d / 2 1 



2 x ttf))- 



(3.13) 



In (S.li) we denote by X_t,JC_j i/ie components of the inverse characteristic 
flow which satisfies $_( o $ t = id; 

*_*(x,0 - [X- t (z,£),*-tM]- 

T7ie associated functions \,v,[i : M + — > R are defined by the following expres- 
sions: 



[D qq a 2 {s) + D pp f3 2 (s) + 2D pq a(s)[3(sj) ds, 
[D qq a 2 {s) + D pp [3 2 (s) + 2D pq a(s)l3(s)^ ds, 



X(t) 



p(t) := 2 / ( D qq a(s)a(s) + D pp f3(s)f3(s) + D pq — (a(s)f3(s)) ) ds, 



where the functions a, (3 : K + — > R are given by: 

a(t) := d- 1 div x (X_ t (z, 0) , P(t) := cT 1 div ? {X- t {x, Q) . 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



Proof. The proof is similar to the calculations given in [Bo|, [Ho]. First note that 
by definition of a and [3, we c an write X- t (x, £) = ot{€)x-\- [3(t)£ . Thus if G is the 
fundamental solution of (3.2), the linear transformation ( 3.1 2| ) guarantees that 
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F has to be a fundamental solution of the following PDE with time dependent 
coefficients 

dtF(t,x,0 = (t)A x + ^(t)A s - ^(t)V a • V S ) F(t,x,0- 

where the functions A, v, /i are calculated depending on the choice of 7. A Fourier 
transform now shows that 

{TF) (t,k,r})=F(t,k,r]):= [ I F{t 7 x, ^) e - l ^ k+ ^ dxd£ 
is a solution of 

d t F(t, k, V ) = - (f Wlfcl 2 + ^WM 2 - • r?)) F(t, fc,7y). 

This equation can easily be integrated and thus, using that F(t, 0, 0) = 1 for all 
t > 0, we obtain 

In (F(t, fc, r?)) = - (A(t)|fc| 2 + v(t)\r,\* - fi(t)(k ■ r?)) . 



After some lengthy and tedious calculations (where one checks that 4\v > n 2 ), 
an inverse Fourier transform gives 

F(t,x,t) = (2ir)- M [ [ e lix k+ ^F(t,k,n)dkdr) 



(47rA) d / 2 J Rd 



( u(t)\x\ 2 +X(t)\j 
\ 4A(t)«/(t) 



LX i J I 4AftWi)-^(t) 



(27r) d (4A(i)K*) - M 2 (*)) rf/2 ' 
which is the desired result. □ 

Note that the quantum mechanical effects in F and consequently in G only 
enter in form of the constants D qqi D pq ~ h 2 , which appear in the auxiliary 
functions A, is, fi. In other words we obtain Green's function for the classical 
Fokker-Planck equation in a square-well potential by setting D qq — D pq = in 
the above expressions. The pointwise positivity of G is a consequence of the 



minimum principle for parabolic equations of Fokker-Planck type [ Ev 



From the above proposition we draw the following consequences (among which 
we obtain the conservativity of the quantum dynamical semigroup, associated 
to the harmonically confined QFP equation). 



Corollary 3.1. Let 7 > and assume (Al)-(AS ). Then for every initial 
condition Wo G L p (M. d x R|) with 1 < p < 00, there exists a unique classical 
solutionwe C(R^;LP(RixRf)) n C 1 (M+ ; C°° (Mf X R|) ) unthM(t) = M = 1. 
Moreover if Wq, no are non-negative a.e., so are w(t, ■), n(t, ■), for all t S IR + . 

Proof. For convenience we use the notation: y := (x, £) as well as yo '■= (xo, £o)- 
With the following linear change of variables 

G(t,$ t {y),y ) = e 2d ^F(t,y-y ) G C 1 (M+, S(R 2d )), 
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we can express our solution in the form 



w(t,$ t {y)) = e 2d ^(w * F(t,y)). 

A straightforward computation now shows that the Jacobian determinant of the 
mapping is given by 



det 



exp(— 2djt). 



With these preparations and using Young's inequality [LiLo], we obtain 

II Mt) \\ p = \\wo* F(t) \\ p < || wo \\ p || F(t) || a < oo, 
since for each fixed t G R + : F(t, •) € S(R 2d ) C L 1 (M 2d ). More precisely we have 



\Ht) 



\F(t,v)\* dy=p-\ 



for all 1 < p < oo. For p = 1 this implies, after a simple calculation, that 
M(t) = Mo = 1, by (A2). Since G(t, •) is pointwise positive, we clearly obtain 
that if wo, resp. no are a.e. non-negative, so are w(t, •), n(t, •). □ 

Remark. If the initial condition wq is the Wigner transform of a pure quantum 
state ipo-, it is well known [Hu|, |LiPa| that wq > pointwise, if and only if ipQ is 
a Gaussian. A similar characterization for mixed states has not been found yet. 

From the above result it is easy to deduce the conservativity (i.e. conservation of 
hcrmiticity, positivity and normalization of the density matrix p) of the quantum 
dynamical semigroup corresponding to ([3.2|) , using the inverse Wigner transform 



p(x,y,t) = (2?r) 



-d 



x + y 



(3.18) 



which is defined in the sense of the usual L 2 -Fourier transform 



| LiLo | . 



3.2 Stationary states. 

As in the classical case, we expect that the competing effects of the confining 
potential plus the positive friction and the dissipating behavior of the operator 
result in a thermal equilibrium state as t — > oo. 

Definition 3.1. A (thermal) equilibrium state is a steady state, i.e. stationary 
solution w € L 2 (R d x K|) of the QFP equation (2.1 ), such that 



J(x)= [ {i-D qq V x ) w{x,£)d£ = 0. (3.19) 

In principle it could be possible that there exist stationary solutions of the QFP 
equation, which are not equilibrium states, however the following proposition 
and corollary show that this is not the case. 

Furthermore, we shall see that the physical intuitive assumption, that the mass 
of the initial state is equal to the mass of the steady state, guarantees its unique- 
ness. 
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Proposition 3.2. Let 7 > 0, ujq > and assume (Al), (Ac ), then the unique 



solution Woo of the stationary QFP equation, satisfying 

Moo '■= [ [ Woo{x,0 dxd£ = l, (3.20) 

JR d JR d 

is given by the following non-isotropic Gaussian function 

Woo(x,£) = expf- — Qhuj 2 \x\ 2 + 2Qi 2 uj x ■ £ + Q 22 \^ 



(2nyVQ V Q 

(3.21) 

such that Woo € S(R d x K^), Q := Q11Q22 — Q\ 2 and 

Qu~D pp + uj%D qq , (3.22) 
g 12 := 2wo7-D M > ( 3 - 23 ) 
Q22 := D pp + LolD qq + A 1 (D pq + 1 D qq ). (3.24) 

Proof. The proof is much simpler in Fourier-space where the steady state Woo 
is explicitly given by 

tMM) = ex P (-£s> (ifcf + ugM 9 )- % |fc| 2 - 
^[(^) N 2 + ^M 2 + ^] 

First it is straightforward to check that this is indeed a solution of the Fourier 
transformed stationary QFP equation 



Lu 2 r) • VfcWoo + (2777 - k) ■ V^Woo = - (D pp \rj\ 2 + D qq \k\ 2 + 2D pq k ■ 77) 



and then that it satisfies the Fourier transformed mass normalization condition 
// Woodxd^ = 1, i.e. 

^00(0,0) = !. (3.25) 

Thus, it remains to prove the uniqueness of the steady state. To do so we 
consider the corresponding characteristic system 

k = u) 2 n, k(s = 0) = fc (3.26) 

77 = 2777 - fc, 77(5 = 0) =770. (3.27) 

Having in mind that, by assumption, 7 > we check that the real parts of 
the eigenvalues Aj.,2 := 7± \/l 2 — w 2 , °f this system are positive, which, by the 
standard theory of ODE's, implies that (k, 77) = (0, 0) is a source of the charac- 



teristic flow (a similar argument holds if 7 = ujq). Thus the condition (3.25) is 



necessary and sufficient to guarantee the uniqueness of the solution Woo ■ d 

From the above proposition we draw the following consequences: 
Corollary 3.2. Under the same assumptions as above we have: 

1. The state Woo is also an equilibrium state, i.e. Joo(x) = for all x € M. d . 
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2. If either loq = (free motion case) or 7 = (frictionless case) or both are 
equal to zero, no nontrivial ^-steady state exists. 

Proof. It is a lengthy but straightforward calculation to show that the current- 
density associated to vanishes identically. 

In the limiting cases 7 = or loq = 0, we consider the characteristic curves given 
by the ODE system (3.26), (3.27). Along these curves the steady state varies 
according to 



-^-M>oo(fc s ,77 s ) = -1? s (fc ,?7o) Wooiks^s), 

as 



(3.28) 



where n s = r/ s (fco , 770) , k s — k s (ko,r]o) denote the (vector valued) characteristic 
curves starting for s = at the point (fco, By ^(fco, 770), we mean 

i?a(fco,fto) := D PP \Vs(ko,Vo)\ 2 + D qq \k s (k ,ri )\ 2 + 2D pq (k s ■ r? s )(fco, 770). 



We integrate equation (3.28) and obtain the Fourier transformed steady state 
parametrized by s € K 



Mko,r)o) dr). (3.29) 



Woo(ks(ko,ri ),r] s (k ,7j )) = M>oo(fco, r?o) exp 



In case the potential vanishes, i.e. oj = (and 7 > 0), we have for all s G 
that fc s (fco,??o) — (free motion) and 



770 - fc T 



if 7 = 0, 



^0 exp(2 7 r) - if 7 >0. 



Now let loq = 7 = and set fc = in equation ( |3.2E ) . We obtain 
Woo(0,77o) =Woo(0,?7o)exp(-£ipp|77o| 2 s), V 770 €E d ,seR, 

which implies 

woo (0, ryo) =0, Vr7 eM d , 



in contradiction to (3.25). The same type of argument holds in the case luq = 0, 
7 > 0, if we set 770 = fco/27 in equation (3.29). 

In the frictionless case, i.e. 7 = 0, ujo > 0, the characteristics are circles in 
the (A;, 77) plane. Thus for s = 2tt/ujq we have (k s (ko,r]o),r] s (ko,rio)) = (kp , 770). 
However if we integrate over one period, i.e. setting s = in ( |3.29 ), we 

obtain for all fco, 770 € M. d 



Woo(k ,ri ) = 77)00 (fc , 770) exp 



;{D vp + ulD qq ) (|fco| 2 +w 2 M 



which clearly implies, since D pp ,D qq > 0, that Woo identically vanishes. □ 

The above corollary in particular shows that, as expected, the presence of a 
positive friction together with a confining potential are crucial to guarantee the 
existence of an admissible, i.e. with finite mass, equilibrium state. 
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Note that although the solution w of ( |2.5| ) in general will not be nonnegative, 
the steady state w M is nevertheless a pointwise positive function, because of 
Assumption ( A2). The associated density matrix is a particular example of 
a mixed quantum state with positive Wigner transform. It is explicitely given by 



poo(x,y) 



7^0 



(3.30) 



Observe that the equilibrium density n^, which is obtained from the diagonal 
(i.e. x = y) of the steady state density matrix Poo ; is a real valued Gaussian 
function. 



Remark. As a special limit case, we observe that in the classical limit, i.e. 

0, the equilibrium state simplifies to the well known station- 



D. 



<i<i 



D 



ary solution of the classical kinetic Fokker-Planck equation (2.2C) with a 
harmonic oscillator potential (c.f. |Ri|), i.e. 



W. 



i(*>0 



w 7 



(IJI 2 +"oM 2 ) 



(2TT) d D p 



(3.31) 



Note that, in contrast to the classical case, the quantum steady state is not a 
function of the classical energy 

H{x,i):=\\tf + \ul\x\\ (3.32) 

For the classical kinetic FP equation, the fact that = w^(H), implies that 
the transport and the classical FP scattering operator vanish independently 
when applied to w^. In our quantum mechanical framework this is no longer 
true, since there the steady state Woo results from a cancellation of the transport 
operator and the scattering term L q . 



4 Long time behavior 

We are now in the position to describe the long time behavior of the linear QFP 
equation in the cases were both friction and an external potential of harmonic 
oscillator type are present, and in cases where one is missing. We will start with 
the latter situation. 



4.1 Dispersion estimates in the unconfined case 

We want to address the unconfined or dispersive cases, i.e. either 7 = or 
c^o = (or both). By comparison with the classical FP equation, we expect 
that the particles escape to infinity and thus that the macroscopic density n 
decays to as t — > 00. More precisely we have the following theorem. 



Theorem 4.1. Let be either 7 = orujQ = 0, or both, and assume (Al). (Ac ): 
1. If Wo £ L 1 H L 2 , then the solution w of the QFP equation satisfies 

II w(t) \\ p < C p R w {ty^ || wo || x , l<p<(X3, (4.1) 
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where C p is a constant independent of Wq, p 1 + q 1 = 1 



R w (t) 



,-4 7 t 



(4\(t)v(t)-^ 2 (t)), 



which is a pointwise positive function with R w (t) — > oo as t — > oo . 
,9. Consequently we have for the corresponding density, that 

II n(t) \\ p < CpRnity^ || wq Hi , 1<p<oo, 
where again the rate R n (t) — > oo as t — > oo. Explicitly 

R n (t) := 2 (A(i)5 2 (t) + fi(t)a(t)0(t) + ^(t)/3 2 (t)) , 

using 

a(t) := cT 1 div xo (X t {x ,&)), 0(f) := d~ x div €o (X t (z 0) £o)). 



(4.2) 

(4.3) 
(4.4) 



Proof. The estimate in claim no . 1 is trivia l in t he case p = oo from the formula 
of the fundamental solution in (3.12) and (2.21). The estimate for p = 1 is due 
to the following property of the Green's function 



// 



G(t,x,£,x ,£ )dx Q d£ = 1 



and a trivial estimate over (2.21). The estimates for 1 < p < oo are then 
obtained by interpolation. The function R w (t) can be computed explicitly in 
each of the three cases and is given by: 



1. If luq > 7 = 0, then with <p = ip(t) — 2 two 
R w (t) 



~ \ D 2 qq + -^D 2 pp ) (^ 2 4 2< us v - 2) 
2 



2 ^pg (C0S<£- 1) 




.Dpp -D g g (cp 2 - 2 cos <^ + 2) 



2 ^pp^qq 
Zuj 



2. If 7 > w = 0, then with % = x(f) = e -2 ' 7 , 
1 



R w (t) 



|7 lPpp + 4 7 Z?pp^)(i 7 (l- X 2 )-(l-x) 2 ) 



1 



7 



2 p<?" /V/ ' ^ ^qq 



D qq (1 - X 2 



3. If 7 = wo = 0, then 
R w (t) - 



4D 2 pq t 2 + AD pp D qq t 2 



- D 2 t 4 

3 ^pp 1 



In all cases, R w diverges as t — > oo and thus claim no. 1 is proved. 

To prove claim no. 2 we first compute the integral of the fundamental solution 
w.r.t. £, which gives 



Gfrx^xo^dZ =R- d / 2 (t)N[ X X ^ o) 

\ y/Rn(t) 



(4.5) 
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where X t is defined as in lemma 3.1 and 

T 2 



Af(v):= (2^)- d / 2 exp^^. 
We can then deduce a formula for the evolution of the macroscopic density, 



using (3.12), (2.21) and the above computation (|4.5|), to obtain 



J J«. d xK. d \ \JRn\t) I 



The dispersion estimates are then straightforward from this expression. On the 
other hand, the function R n (t) again can be computed explicitly: 

1. If uq > 7 — 0, then with <p — (p(t) — 2 two, 

2 R n (t) — 2ui D pq (1 — cos (pj + D pp (ip — sin^) + u 2 D qq (<p + shxip) . 

2. If 7 > uj = 0, then with \ = X(*) = e -2 * 7 , 

8j 3 R n (t) = D pp (A X + 4 7 i - 3 - x 2 ) + l&ftD m + 8jD pq ( X + 2 7 t - 1) . 

3. If 7 — ujq = 0, then 

Rn (t) = 2 D g9 < + 2 £> p , t 2 + | Dpp t 3 . 

In all three cases R n diverges as t — > oo and thus claim no. 2 is also proved. □ 

The behavior at t = +oo of the rate functions can be obtained directly from 
the explicit computations in the previous result, depending on the diffusion 
constants. In the following corollary we will consider only the physical most 
important case, namely D pp > 0. 



Corollary 4.1. Assume D pp > and (Al), ^Aoj ), then we have, as t — > +oo: 

1. If luq > 7 = 0, then R w (t) = 0(t 2 ) and R n {t) = 0(t). 

2. If 7 > ujq = 0, then R w (t) = 0(t) and R n (t) = 0(t). 

3. If 7 = ojq = 0, then R w (t) = 0{t 4 ) and R n {t) = 0(t 3 ). 



Remark. The computation of n can also be used to measure the convergence 
rates towards the steady state macroscopic density in the confined case. 
This computation however is quite involved and we leave the details to the 
reader, since we will take a more elegant approach in the next subsection. 

Nevertheless the advantage of such an explicit calculation would be that one 
can prove exponential decay of the solution, even for the classical FP equation 
(this can also be seen by comparison with spectral theoretical approaches p^|), 
whereas the entropy method used in section 4.2 only gives a suboptimal rate in 



this case [DeVi 
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4.2 Exponential decay towards equilibrium 



We now assume the presence of friction and of the confining potential. Like 
in the classical case, we expect exponential decay of the solution towards the 
equilibrium state, which is usually proved using spectral theory. 

However in this work we shall follow a different approach, the so-called entropy- 
entropy- dissipation method for classical FP type equations (see [AMTU] and 
references therein). As we shall see it can be successfully applied in our quantum 
mechanical context too (for an overview of the classical applications, see [ MaVfl ) . 

First let us rewrite the drift-vector P given by (3.6) in the following way 

P(x,0=D(VA(x,^)+F(x,0), (4.6) 

where the gradient is taken w.r.t. (x, £) £ M. 2d and A is defined as the normalized 
potential appearing in the expression of the equilibrium state, more precisly 



(4.7) 



Consequently the QFP equation ( |3.2| ) takes the standard form of a non-sym- 
metric drift- diffusion equation 



d t w = div [D(Vw + w(VA + F))] , t £ M+ 
w(t = 0,x,f) = w (x,£). 

The vector field F is explicitly given by 
F(x,£) :=D- 1 P(x,0-VA(x,0, 



(4.8) 
(4.9) 

(4.10) 



where D 1 is the i nve rse of the diffusion matrix, A is defined by (4/7) and P is 
the drift given in (3.6). 

The rewritten QFP e qua tion (for harmonic oscillator potentials) with equilib- 
rium state giv en by (4.7 ) can now be identified as a special case of FP type 
equations, see [AMTUJ. Therefore we introduce the concept of relative entropies 
in the same sense as in the quoted work: 

Definition 4.1. Let A be either R orR+. Let ip £ C{A) n C A {A) satisfy 

(p(l) = 0, if" > with if" ^ and ((p'") 2 < -<p"(p IV on A. (4.11) 

Assume moreover that f £ L 1 (IR™), g £ L+(M. n ) with J R „ f dx — L„ g dx = 1 
and f jg £ A g{dx) a.e. Then 



e v {f I 5) 



( ^ ) .'/('/•'•') 



(4.12) 



is an admissible relative entropy of f w.r.t. g and generating function ip. 

Two relative entropies which are frequently used are the logarithmic relative 
entropy, associated to the generating function 



(pi(ct) := alna — a + 1, a £ R + . 
and the quadratic relative entropy with generator 
ip 2 (a) := k(a - l) 2 , a £ E, k> 0. 



(4.13) 
(4.14) 
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Indeed it is known, see [AMTU], that (up to a positive multiplicative constant) 
for every admissible entropy generator ip, there exist constants K\, Kq, > such 
that for all a £ R + , it holds that 



Kiipt(a) < <p(a) < K 2 ip2{a). 



(4.15) 



Remark. The above dehned entropies should not be confused with the quantum 
mechanical von Neumann entropy S := — Tr(i?lni?), where R is the density 
operator of the particle ensemble |Th[ (also note that in contrast to the physical 
convention the minus sign in front of (4.12) is dropped). 

There is however a notion of entropy for quantum states, called Wehrl entropy 
[ We ] , which is closely related to the logarithmic entropy defined above (using 
the pointwise nonnegative Husirai transform of w, see, e.g., [ LiPaj ) and which 
can be interpreted as a measure of coherence and localization of quantum states. 



For details, sec [AnHa], [GnZy|, [SIZy] and the references given therein. 

Since the solution of the QFP equation in general is not pointwise positive, it 
seems that wc need to look for an admissible entropy on all of K. This would 



imply, see [AMTU], that the the only admissible entropy useful for our purpose 
is the quadratic one (4.14). We circumvent this shortcoming by decomposing 
the initial condition wq into its negative and positive parts. More precisely we 
write 



w Q (x, = w£(x, f) - w (x, £) 



(4.16) 



where now Wq , w are both non-negative functions with mass Mq , M respec- 



tively. Clearly condition (A2) implies 



l = M n 4 



Afn 



Having in mind that the QFP equation is linear, we denote by 



w ± (t,x,£) :-- 



Wq(xo,€o) G(t,x, £,x ,£o) dx d£ 



(4.17) 



(4.18) 



the mild solution of (4.8) corresponding to iujj\ Corollary 2.1 implies that 
(t, •) > a.e. and we shall now apply the entropy-entropy-dissipation method 
to each of the two functions w + , w~ . Note that the steady state associated to 
w ± (t, •) is given by = MqW^. 

In order to use the results established in [AMTU], we first need to check that 
the following property is fulfilled. 

Lemma 4.1. Assume 7 > 0, loq > as well as ( |A y := (x, £) € 
Then it holds that 



p2rf 



divy (DFwoo) — on 



],2d 



(4.19) 



where F is defined by ( (.1C ) 



Proof. To prove the claim, we note, using equation (4.7), that condition (4.19) 
is equivalent to 

divy (DVwoo + Pwoo) = 0, 



which is the stationary QFP equation for potentials of the form (A3). Thus the 
claim is true by definition of A, F and Woq. □ 
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With these preparations we can now state the main result of this section. 



Theorem 4.2. Assume 7 > 0, loq > and assume (Al)-(Ah). Let the initial 
data wo be such that Wq(x,^) = Wq(x,^) — Wq(x,^) a.e., w G L 1 n L 2 and 
assume 



e^O^ I M ± w oc ) < 



Then there exists a k > such that 



e ip {w ± {t,-) I MqW, 



< e 



-2/tt, 



e^(u^ I M ± it? oo ), t > 0. 



(A4) 



(4.20) 



As a consequence, the solution of (S.i) converges exponentially towards the equi- 
librium state. More precisely it holds that: 



1. If e Vl (u> I Woo) < 00, with ipi defined in (4-L; ), then 



II 



w c 



li <Ce 



C e R + , t> 0. 



(4.21) 



2. More generally, if e ¥ , 2 (wg | w^,) < 00, with ip 2 defined in then 

II w(t,-) - Woo \\ p < Ce- Kt CtR+,t>0, (4.22) 
with 1 < p < 2. 



Proof. The proof of the first claim (4.2C ) is a consequence of the following convex 
Sobolev inequality AMTU (in which we use the notation dpoo = p^ dx) 



tp ( — I dpoz < ^— 

Po ■ 



9 



P 

Poo 



p 

Poo 



dpoo, 



which holds for every admissible entropy generator ip (c.f. definition 4.1), every 
function p € L 1 (R 2d ), with M = and Poo := exp(-A(x, £)) € (M 2d ), 
with uniformly convex potential A, i.e. 



d 2 A 



where K\ > and l2d denotes the identity matrix in 
we can now estimate (recall that = M ± w 00 ) 



(4.23) 

i 2d . Using this inequality 



<P 



^\dwt< 



1 



25m 



j.2,; 



dwt 



where 6 > is the smallest eigenvalue of D given by 
5 := 



D„ 



D ) 2 



AD 2 

pq 



(4.24) 



sin ce by defi nition of 5, it holds that Sl2d < D. Provided lemma 4.1, the results 
of I AMTU I imply the exponential decay of the relative entropy with a rate 
k = Ski and thus ( 4.20| ) is proved. 



Long time behavior of QFP 



18 



Using the well known Csiszdr-Kullback inequality |AMTU | , [ Cs | , | Ku | , we further 
obtain 



II w ± (t,-) 



M^uUL < Ke 



L "'ooll! 

Since by (A2) we have 1 = M t 



~ K \ KeM+,t>0. 
AIq , this implies 



\w(t, •) - Woc || 1 = \\w+(t, .)-«;-(*,.)- (Af+ - Mo)w oa \\ 1 

< \\w + (t,-) -Mo+Woolli + ||t«-(t,-) - AfJ-lOooHi 

< Ce~ K \ CeR + . 



For entropy generators of the form ( 4.14 ) , we get from ( 4.2Cj ) 

|| w ± (t, •) - M ± u; oo || 2>5 < e- Kt || w$ - M± Woo \\ 2S Vt>0, 

where 



| 2 s is the L -norm with weight S = l/u>oo, i.e. 



/II 



2,5 



R d xR d 



As above, the use of the triangle inequality allows us to conclude that exponen- 
tial decay (with rate n) of w holds in this weighted L 2 -norm. Since l/woo > 1 
poin twise, we have |j / || 2 < || / || 2 g (here we have used the Lindblad condition 



( A.1) to show that the constant in front of l/woo is indeed greater than 1). 
This implies exponential decay in the usual L 2 -norm and thus, by interpolation, 
we obtain the exponential convergence of w towards the steady state in all 
Z/ p -norms with 1 < p < 2. □ 



It should be noted, that the main drawback of the above theorem is the fact 
that it only holds for initial data with positive and negative part bounded in 
relative entropy (A4). Furthermore one should have in mind that the above 
theorem fails, if 7 = or loq = 0, since then A is no longer uniformly convex. 
However we have already seen that in these cases no finite-mass steady state, 
different from zero, exists. 

The precise value of the rate k can be obtained by the following transformation. 
Define a new potential A by 



:=A IVD (ar.0 



(4.25) 



where VD is the square root of D in the sense of positive definite matrices. 
Then, it holds that 



d 2 A 



(4.26) 



which gives the optimal rate k as the smallest eigenvalue of the Hessian of A. 



For details, see ]AMTU| or (RJ. 

Remark. In contrast to our result, the convergence to the equilibrium state for 
classical kinetic FP equations ( 2.20| ) is, by similar entropy methods, obtained 
with a suboptimal rate of order 0(i~°°), see [DcVi|. This is due to the fact, 
that for the classical FP equation, in addition to the thermal equilibrium state 
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( |3.31 ), there exist other states with zero entropy dissipation. Indeed explicit 
calculations show that k depends on H and converges to as k — > 0. 
We hope to extend the above results to more general, maybe nonlinear, poten- 
tials V in forthcoming works. Although quite simple, the presented harmonic 
oscillator case is expected to be an essential prerequisite in this analysis. 
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